An immunoinformatics approach to study the epitopes of SARS-CoV-2 helicase, Nsp13.

Introduction and objective: Vaccines are administered worldwide to control on-going coronavirus disease-19 (COVID-19) pandemic caused by SARS-CoV-2. Vaccine efficacy is largely contributed by the epitopes present on the viral proteins and their alteration might help emerging variants to escape host immune surveillance. Therefore, this study was designed to study SARS-CoV-2 Nsp13 protein, its epitopes and evolution. Methods: Clustal Omega was used to identify mutations in Nsp13 protein. Secondary structure and disorder score was predicted by CFSSP and PONDR-VSL2 webservers. Protein stability was predicted by DynaMut webserver. B cell epitopes were predicted by IEDB DiscoTope 2.0 tools and their 3D structures were represented by discovery studio. Antigenicity and allergenicity of epitopes were predicted by Vaxijen2.0 and AllergenFPv.1.0. Physiochemical properties of epitopes were predicted by Toxinpred, HLP webserver tool. Results: Our data revealed 182 mutations in Nsp13 among Indian SARS-CoV-2 isolates, which were characterised by secondary structure and per-residue disorderness, stability and dynamicity predictions. To correlate the functional impact of these mutations, we characterised the most prominent B cell and T cell epitopes contributed by Nsp13. Our data revealed twenty-one epitopes, which exhibited antigenicity, stability and interactions with MHC class-I and class-II molecules. Subsequently, the physiochemical properties of these epitopes were analysed. Furthermore, eighteen mutations reside in these Nsp13 epitopes. Conclusions: We report appearance of eighteen mutations in the predicted twenty-one epitopes of Nsp13. Among these, at least seven epitopes closely matches with the functionally validated epitopes. Altogether, our study shows the pattern of evolution of Nsp13 epitopes and their probable implications.

Introducción y objetivo: Las vacunas se administran a nivel mundial para controlar la pandemia en curso de la enfermedad por coronavirus de 2019 (COVID-19) causada por SARS-CoV-2. A la eficacia de la vacuna contribuyen ampliamente los epítopes presentes en las proteínas virales, y su alteración puede contribuir a que las variantes emergentes se escapen de la vigilancia inmunológica del huésped. Por tanto, este estudio fue diseñado para estudiar la proteína Nsp13 de SARS-CoV-2, sus epítopes y su evolución. Métodos: Se utilizó Clustal Omega para identificar las mutaciones de la proteína Nsp13. La estructura secundaria y la tasa de desorden se predijeron mediante los servidores web CFSSP y PONDR-VSL2. La estabilidad de la proteína fue predicha mediante el servidor web DynaMut. Los epítopes de las células B fueron predichos mediante las herramientas DiscoTope 2.0 de IEDB, y sus estructuras en 3D fueron representadas mediante Discovery Studio.La antigenicidad y alergenicidad de los epítopes fueron predichas mediante Vaxijen2.0 y AlergenFPv.1.0. Las propiedades fisioquímicas de los epítopes fueron predichas mediante Toxinpred, la herramienta del servidor web HLP. Resultados: Nuestros datos revelaron 182 mutaciones en Nsp13 entre los aislados indios de SARS-CoV-2, que fueron caracterizadas mediante las predicciones de la estructura secundaria y la capacidad de desorden por residuo, la estabilidad y la dinamicidad. Para correlacionar el impacto funcional de estas mutaciones, caracterizamos los epítopes más prominentes de las células B y las células T a los que contribuyó Nsp13. Nuestros datos revelaron veintiún epítopes, que exhibieron antigenicidad, estabilidad e interacciones con las moléculas MHC de clase I y clase II. Seguidamente se analizaron las propiedades fisioquímicas de estos epítopes. Además, en estos epítopes de Nsp13 residen ocho mutaciones. Conclusiones: Reportamos el aspecto de ocho mutaciones en los veintiún epítopes de Nsp13 predichos. Entre estos, al menos siete epítopes concuerdan estrechamente con los epítopes funcionalmente validados. En su conjunto, nuestro estudio refleja el patrón evolutivo de los epítopes de Nsp13 y sus implicaciones probables.

Non-uniform aspects of the SARS-CoV-2 intraspecies evolution reopen question of its origin.

Several hypotheses have been presented on the origin of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) from its identification as the agent causing the current coronavirus disease 19 (COVID-19) pandemic. So far, no solid evidence has been found to support any hypothesis on the origin of this virus, and the issue continue to resurface over and over again. Here we have unfolded a pattern of distribution of several mutations in the SARS-CoV-2 proteins in 24 geo-locations across different continents. The results showed an evenly uneven distribution of the unique protein variants, distinct mutations, unique frequency of common conserved residues, and mutational residues across these 24 geo-locations. Furthermore, ample mutations were identified in the evolutionarily conserved invariant regions in the SARS-CoV-2 proteins across almost all geo-locations studied. This pattern of mutations potentially breaches the law of evolutionary conserved functional units of the beta-coronavirus genus. These mutations may lead to several novel SARS-CoV-2 variants with a high degree of transmissibility and virulence. A thorough investigation on the origin and characteristics of SARS-CoV-2 needs to be conducted in the interest of science and for the preparation of meeting the challenges of potential future pandemics.

An issue of concern: unique truncated ORF8 protein variants of SARS-CoV-2.

Open reading frame 8 (ORF8) shows one of the highest levels of variability among accessory proteins in Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), the causative agent of Coronavirus Disease 2019 (COVID-19). It was previously reported that the ORF8 protein inhibits the presentation of viral antigens by the major histocompatibility complex class I (MHC-I), which interacts with host factors involved in pulmonary inflammation. The ORF8 protein assists SARS-CoV-2 in evading immunity and plays a role in SARS-CoV-2 replication. Among many contributing mutations, Q27STOP, a mutation in the ORF8 protein, defines the B.1.1.7 lineage of SARS-CoV-2, engendering the second wave of COVID-19. In the present study, 47 unique truncated ORF8 proteins (T-ORF8) with the Q27STOP mutations were identified among 49,055 available B.1.1.7 SARS-CoV-2 sequences. The results show that only one of the 47 T-ORF8 variants spread to over 57 geo-locations in North America, and other continents, which include Africa, Asia, Europe and South America. Based on various quantitative features, such as amino acid homology, polar/non-polar sequence homology, Shannon entropy conservation, and other physicochemical properties of all specific 47 T-ORF8 protein variants, nine possible T-ORF8 unique variants were defined. The question as to whether T-ORF8 variants function similarly to the wild type ORF8 is yet to be investigated. A positive response to the question could exacerbate future COVID-19 waves, necessitating severe containment measures.

Emergence of unique SARS-CoV-2 ORF10 variants and their impact on protein structure and function.

The devastating impact of the ongoing coronavirus disease 2019 (COVID-19) on public health, caused by the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has made targeting the COVID-19 pandemic a top priority in medical research and pharmaceutical development. Surveillance of SARS-CoV-2 mutations is essential for the comprehension of SARS-CoV-2 variant diversity and their impact on virulence and pathogenicity. The SARS-CoV-2 open reading frame 10 (ORF10) protein interacts with multiple human proteins CUL2, ELOB, ELOC, MAP7D1, PPT1, RBX1, THTPA, TIMM8B, and ZYG11B expressed in lung tissue. Mutations and co-occurring mutations in the emerging SARS-CoV-2 ORF10 variants are expected to impact the severity of the virus and its associated consequences. In this article, we highlight 128 single mutations and 35 co-occurring mutations in the unique SARS-CoV-2 ORF10 variants. The possible predicted effects of these mutations and co-occurring mutations on the secondary structure of ORF10 variants and host protein interactomes are presented. The findings highlight the possible effects of mutations and co-occurring mutations on the emerging 140 ORF10 unique variants from secondary structure and intrinsic protein disorder perspectives.

Emerging genetic diversity of SARS-CoV-2 RNA dependent RNA polymerase (RdRp) alters its B-cell epitopes.

The RNA dependent RNA polymerase (RdRp) plays crucial role in virus life cycle by replicating the viral genome. The SARS-CoV-2 is an RNA virus that rapidly spread worldwide and acquired mutations. This study was carried out to identify mutations in RdRp as the SARS-CoV-2 spread in India. We compared 50217 RdRp sequences reported from India with the first reported RdRp sequence from Wuhan, China to identify 223 mutations acquired among Indian isolates. Our protein modelling study revealed that several mutants can potentially alter stability and flexibility of RdRp. We predicted the potential B cell epitopes contributed by RdRp and identified thirty-six linear continuous and twenty-five discontinuous epitopes. Among 223 RdRp mutants, 44% of them localises in the B cell epitopes region. Altogether, this study highlights the need to identify and characterize the variations in RdRp to understand the impact of these mutations on SARS-CoV-2.

Implications derived from S-protein variants of SARS-CoV-2 from six continents.

The spike (S) protein is a critical determinant of the infectivity and antigenicity of SARS-CoV-2. Several mutations in the S protein of SARS-CoV-2 have already been detected, and their effect in immune system evasion and enhanced transmission as a cause of increased morbidity and mortality are being investigated. From pathogenic and epidemiological perspectives, S proteins are of prime interest to researchers. This study focused on the unique variants of S proteins from six continents: Asia, Africa, Europe, Oceania, South America, and North America. In comparison to the other five continents, Africa had the highest percentage of unique S proteins (29.1%). The phylogenetic relationship implies that unique S proteins from North America are significantly different from those of the other five continents. They are most likely to spread to the other geographic locations through international travel or naturally by emerging mutations. It is suggested that restriction of international travel should be considered, and massive vaccination as an utmost measure to combat the spread of the COVID-19 pandemic. It is also further suggested that the efficacy of existing vaccines and future vaccine development must be reviewed with careful scrutiny, and if needed, further re-engineered based on requirements dictated by new emerging S protein variants.

Periodically aperiodic pattern of SARS-CoV-2 mutations underpins the uncertainty of its origin and evolution.

Various lineages of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) have contributed to prolongation of the Coronavirus Disease 2019 (COVID-19) pandemic. Several non-synonymous mutations in SARS-CoV-2 proteins have generated multiple SARS-CoV-2 variants. In our previous report, we have shown that an evenly uneven distribution of unique protein variants of SARS-CoV-2 is geo-location or demography-specific. However, the correlation between the demographic transmutability of the SARS-CoV-2 infection and mutations in various proteins remains unknown due to hidden symmetry/asymmetry in the occurrence of mutations. This study tracked how these mutations are emerging in SARS-CoV-2 proteins in six model countries and globally. In a geo-location, considering the mutations having a frequency of detection of at least 500 in each SARS-CoV-2 protein, we studied the country-wise percentage of invariant residues. Our data revealed that since October 2020, highly frequent mutations in SARS-CoV-2 have been observed mostly in the Open Reading Frame (ORF) 7b and ORF8, worldwide. No such highly frequent mutations in any of the SARS-CoV-2 proteins were found in the UK, India, and Brazil, which does not correlate with the degree of transmissibility of the virus in India and Brazil. However, we have found a signature that SARS-CoV-2 proteins were evolving at a higher rate, and considering global data, mutations are detected in the majority of the available amino acid locations. Fractal analysis of each protein's normalized factor time series showed a periodically aperiodic emergence of dominant variants for SARS-CoV-2 protein mutations across different countries. It was noticed that certain high-frequency variants have emerged in the last couple of months, and thus the emerging SARS-CoV-2 strains are expected to contain prevalent mutations in the ORF3a, membrane, and ORF8 proteins. In contrast to other beta-coronaviruses, SARS-CoV-2 variants have rapidly emerged based on demographically dependent mutations. Characterization of the periodically aperiodic nature of the demographic spread of SARS-CoV-2 variants in various countries can contribute to the identification of the origin of SARS-CoV-2.

Autoimmunity roots of the thrombotic events after COVID-19 vaccination.

Although vaccination represents the most promising way to stop or contain the coronavirus disease 2019 (COVID-19) pandemic and safety and effectiveness of available vaccines were proven, a small number of individuals who received anti-SARS-CoV-2 vaccines developed a prothrombotic syndrome. Vaccine-induced immune thrombotic thrombocytopenia (VITT) can be triggered by the adenoviral vector-based vaccine, whereas lipid nanoparticle-mRNA-based vaccines can induce rare cases of deep vein thrombosis (DVT). Although the main pathogenic mechanisms behind this rare phenomenon have not yet been identified, both host and vaccine factors might be involved, with pathology at least in part being related to the vaccine-triggered autoimmune reaction. In this review, we are considering some aspects related to pathogenesis, major risk factors, as well as peculiarities of diagnosis and treatment of this rare condition.

Urgent Need for Field Surveys of Coronaviruses in Southeast Asia to Understand the SARS-CoV-2 Phylogeny and Risk Assessment for Future Outbreaks.

Phylogenetic analysis of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is focused on a single isolate of bat coronaviruses (bat CoVs) which does not adequately represent genetically related coronaviruses (CoVs) [...].

COVID-19 Vaccines and Thrombosis-Roadblock or Dead-End Street?

Two adenovirus-based vaccines, ChAdOx1 nCoV-19 and Ad26.COV2.S, and two mRNA-based vaccines, BNT162b2 and mRNA.1273, have been approved by the European Medicines Agency (EMA), and are invaluable in preventing and reducing the incidence of coronavirus disease-2019 (COVID-19). Recent reports have pointed to thrombosis with associated thrombocytopenia as an adverse effect occurring at a low frequency in some individuals after vaccination. The causes of such events may be related to SARS-CoV-2 spike protein interactions with different C-type lectin receptors, heparan sulfate proteoglycans (HSPGs) and the CD147 receptor, or to different soluble splice variants of the spike protein, adenovirus vector interactions with the CD46 receptor or platelet factor 4 antibodies. Similar findings have been reported for several viral diseases after vaccine administration. In addition, immunological mechanisms elicited by viral vectors related to cellular delivery could play a relevant role in individuals with certain genetic backgrounds. Although rare, the potential COVID-19 vaccine-induced immune thrombotic thrombocytopenia (VITT) requires immediate validation, especially in risk groups, such as the elderly, chronic smokers, and individuals with pre-existing incidences of thrombocytopenia; and if necessary, a reformulation of existing vaccines.

Notable sequence homology of the ORF10 protein introspects the architecture of SARS-CoV-2.

The current Coronavirus Disease 19 (COVID-19) pandemic, caused by Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) shows similar pathology to MERS and SARS-CoV, with a current estimated fatality rate of 1.4%. Open reading frame 10 (ORF10) is a unique SARS-CoV-2 accessory protein, which contains eleven cytotoxic T lymphocyte (CTL) epitopes each of nine amino acids in length. Twenty-two unique SARS-CoV-2 ORF10 variants have been identified based on missense mutations found in sequence databases. Some of these mutations are predicted to decrease the stability of ORF10 in silico physicochemical and structural comparative analyses were carried out on SARS-CoV-2 and Pangolin-CoV ORF10 proteins, which share 97.37% amino acid (aa) homology. Though there is a high degree of ORF10 protein similarity of SARS-CoV-2 and Pangolin-CoV, there are differences of these two ORF10 proteins related to their sub-structure (loop/coil region), solubility, antigenicity and shift from strand to coil at aa position 26 (tyrosine). SARS-CoV-2 ORF10, which is apparently expressed in vivo since reactive T cell clones are found in convalescent patients should be monitored for changes which could correlate with the pathogenesis of COVID-19.

A unique view of SARS-CoV-2 through the lens of ORF8 protein.

Immune evasion is one of the unique characteristics of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) attributed to its ORF8 protein. This protein modulates the adaptive host immunity through down-regulation of MHC-1 (Major Histocompatibility Complex) molecules and innate immune responses by surpassing the host's interferon-mediated antiviral response. To understand the host's immune perspective in reference to the ORF8 protein, a comprehensive study of the ORF8 protein and mutations possessed by it have been performed. Chemical and structural properties of ORF8 proteins from different hosts, such as human, bat, and pangolin, suggest that the ORF8 of SARS-CoV-2 is much closer to ORF8 of Bat RaTG13-CoV than to that of Pangolin-CoV. Eighty-seven mutations across unique variants of ORF8 in SARS-CoV-2 can be grouped into four classes based on their predicted effects (Hussain et al., 2021) [1]. Based on the geo-locations and timescale of sample collection, a possible flow of mutations was built. Furthermore, conclusive flows of amalgamation of mutations were found upon sequence similarity analyses and consideration of the amino acid conservation phylogenies. Therefore, this study seeks to highlight the uniqueness of the rapidly evolving SARS-CoV-2 through the ORF8.

Variations in Orf3a protein of SARS-CoV-2 alter its structure and function.

Severe acquired respiratory syndrome coronavirus 2 (SARS-CoV-2) rapidly spread worldwide and acquired multiple mutations in its genome. Orf3a, an accessory protein encoded by the genome of SARS-CoV-2, plays a significant role in viral infection and pathogenesis. In the present in-silico study, 15,928 sequences of Orf3a reported worldwide were compared to identify variations in this protein. Our analysis revealed the occurrence of mutations at 173 residues of Orf3a protein. Subsequently, protein modelling was performed that revealed twelve mutations which can considerably affect the stability of Orf3a. Among the 12 mutations, three mutations (Y160H, D210Y and S171L) also lead to alterations in secondary structure and protein disorder parameters of the Orf3a protein. Further, we used predictive tools to identify five promising epitopes of B-cells, which resides in the mutated regions of Orf3a. Altogether, our study sheds light on the variations occurring in Orf3a that might contribute to alteration in protein structure and function.

Identification of twenty-five mutations in surface glycoprotein (Spike) of SARS-CoV-2 among Indian isolates and their impact on protein dynamics.

SARS-CoV-2, the causative agent of the COVID-19 pandemic, is an RNA virus that has inherent high rate of mutation. Due to the mutations, the virus evolves at a rapid pace that helps them to survive better inside the host. One of the hotspots of pharmacological interventions is to inhibit binding of virus with the host cells, which is mediated by Spike glycoprotein of SARS-CoV-2 and ACE2 receptors present on the human cells. This study was conducted with an aim to identify and characterise the mutation (s) present in the Spike glycoprotein of the SARS-CoV-2. Towards this, an in silico methodology was used, and the mutations on Spike glycoprotein were identified by comparing the Spike glycoprotein of first reported sequence from Wuhan wet seafood market virus with the available sequences of SARS-CoV-2 from Indian isolates. Our analysis revealed the presence of twenty-five mutations in Spike glycoprotein among Indian SARS-CoV-2 isolates. These mutations spread all over the protein and can be clustered at least into four distinct positions. Further, mutations at eleven positions exhibited alterations in the secondary structure of the polypeptide chain. We also investigated the influence of these mutations on overall protein dynamics and have shown that they affect the dynamic stability of the Spike glycoprotein.

Identification and molecular characterization of mutations in nucleocapsid phosphoprotein of SARS-CoV-2.

SARS-CoV-2 genome encodes four structural proteins that include the spike glycoprotein, membrane protein, envelope protein and nucleocapsid phosphoprotein (N-protein). The N-protein interacts with viral genomic RNA and helps in packaging. As SARS-CoV-2 spread to almost all countries worldwide within 2-3 months, it also acquired mutations in its RNA genome. Therefore, this study was conducted with an aim to identify the variations present in N-protein of SARS-CoV-2. Here, we analysed 4,163 reported sequence of N-protein from United States of America (USA) and compared them with the first reported sequence from Wuhan, China. Our study identified 107 mutations that reside all over the N-protein. Further, we show the high rate of mutations in intrinsically disordered regions (IDRs) of N-protein. Our study show 45% residues of IDR2 harbour mutations. The RNA-binding domain (RBD) and dimerization domain of N-protein also have mutations at key residues. We further measured the effect of these mutations on N-protein stability and dynamicity and our data reveals that multiple mutations can cause considerable alterations. Altogether, our data strongly suggests that N-protein is one of the mutational hotspot proteins of SARS-CoV-2 that is changing rapidly and these mutations can potentially interferes with various aspects of N-protein functions including its interaction with RNA, oligomerization and signalling events.

Possible Transmission Flow of SARS-CoV-2 Based on ACE2 Features.

Angiotensin-converting enzyme 2 (ACE2) is the cellular receptor for the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) that is engendering the severe coronavirus disease 2019 (COVID-19) pandemic. The spike (S) protein receptor-binding domain (RBD) of SARS-CoV-2 binds to the three sub-domains viz. amino acids (aa) 22-42, aa 79-84, and aa 330-393 of ACE2 on human cells to initiate entry. It was reported earlier that the receptor utilization capacity of ACE2 proteins from different species, such as cats, chimpanzees, dogs, and cattle, are different. A comprehensive analysis of ACE2 receptors of nineteen species was carried out in this study, and the findings propose a possible SARS-CoV-2 transmission flow across these nineteen species.

The structural basis of accelerated host cell entry by SARS-CoV-2†.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the pandemic coronavirus disease 2019 (COVID-19) that exhibits an overwhelming contagious capacity over other human coronaviruses (HCoVs). This structural snapshot describes the structural bases underlying the pandemic capacity of SARS-CoV-2 and explains its fast motion over respiratory epithelia that allow its rapid cellular entry. Based on notable viral spike (S) protein features, we propose that the flat sialic acid-binding domain at the N-terminal domain (NTD) of the S1 subunit leads to more effective first contact and interaction with the sialic acid layer over the epithelium, and this, in turn, allows faster viral 'surfing' of the epithelium and receptor scanning by SARS-CoV-2. Angiotensin-converting enzyme 2 (ACE-2) protein on the epithelial surface is the primary entry receptor for SARS-CoV-2, and protein-protein interaction assays demonstrate high-affinity binding of the spike protein (S protein) to ACE-2. To date, no high-frequency mutations were detected at the C-terminal domain of the S1 subunit in the S protein, where the receptor-binding domain (RBD) is located. Tight binding to ACE-2 by a conserved viral RBD suggests the ACE2-RBD interaction is likely optimal. Moreover, the viral S subunit contains a cleavage site for furin and other proteases, which accelerates cell entry by SARS-CoV-2. The model proposed here describes a structural basis for the accelerated host cell entry by SARS-CoV-2 relative to other HCoVs and also discusses emerging hypotheses that are likely to contribute to the development of antiviral strategies to combat the pandemic capacity of SARS-CoV-2.

Identification of novel mutations in the methyltransferase complex (Nsp10-Nsp16) of SARS-CoV-2.

A recent outburst of the pandemic caused by a member of the coronaviridae family identified as SARS-CoV-2. The highly contagious nature of the virus allows it to spread rapidly worldwide and caused severe healthcare and economic distress. So far, no proper line of treatment or vaccines has been available against SARS-CoV-2. Since, the infected people rapidly increased, causing the saturation of healthcare systems with coronavirus disease (COVID-19) patients. As the virus spread to new locations it also acquired various mutations. Here, in this study, we focused on identifying mutations in one of the crucial complex of SARS-CoV-2, the Nsp10-Nsp16 2'-O-methyltransferase complex. This complex plays indispensable role in the post-transcriptional modifications of viral RNA by its capping. We analysed 208 sequences of Nsp10-Nsp16 reported from India and compared with first reported sequence from Wuhan, China. Our analysis revealed a single mutation in Nsp10 and five mutations in Nsp16 protein. We also show that these mutations are leading to alteration in the secondary structure of Nsp10-Nsp16. Further, the protein modelling studies revealed that the mutation of both Nsp10-Nsp16 impacts the protein dynamicity and stability. Altogether, this study provides novel insights into the variations observed in the proteins of SARS-CoV-2 that might have functional consequences.

Identification of novel mutations in RNA-dependent RNA polymerases of SARS-CoV-2 and their implications on its protein structure.

The rapid development of the SARS-CoV-2 mediated COVID-19 pandemic has been the cause of significant health concern, highlighting the immediate need for effective antivirals. SARS-CoV-2 is an RNA virus that has an inherently high mutation rate. These mutations drive viral evolution and genome variability, thereby facilitating viruses to have rapid antigenic shifting to evade host immunity and to develop drug resistance. Viral RNA-dependent RNA polymerases (RdRp) perform viral genome duplication and RNA synthesis. Therefore, we compared the available RdRp sequences of SARS-CoV-2 from Indian isolates and the 'Wuhan wet sea food market virus' sequence to identify, if any, variation between them. Our data revealed the occurrence of seven mutations in Indian isolates of SARS-CoV-2. The secondary structure prediction analysis of these seven mutations shows that three of them cause alteration in the structure of RdRp. Furthermore, we did protein modelling studies to show that these mutations can potentially alter the stability of the RdRp protein. Therefore, we propose that RdRp mutations in Indian SARS-CoV-2 isolates might have functional consequences that can interfere with RdRp targeting pharmacological agents.